Suppression of Spin Frustration due to Orbital Selection 
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In order to clarify a crucial role of orbital degree of freedom in geometrically frustrated systems, 
we investigate both ground- and excited-state properties of the e g -orbital degenerate Hubbard model 
on two kinds of lattices, ladder and zigzag chain, by using numerical techniques. In the ladder, spin 
correlation extends on the whole system, while the zigzag chain is decoupled to a double chain and 
spin excitation is confined in one side of the double chain due to the selection of a specific orbital. 
We envision a kind of self-organization phenomenon that the geometrically frustrated multi-orbital 
system is spontaneously reduced to a one-orbital model to suppress the spin frustration. 

PACS numbers: 75.10.-b, 71.10.Fd, 75.30.Et, 75.40.Mg 
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Recently strongly correlated electron systems with ge- 
ometrical frustration have attracted much attention in 
the research field of condensed-matter physics, since a 
subtle balance among competing interactions leads to a 
variety of cooperative phenomena such as exotic super- 
conductivity and novel magnetism. In particular, the 
discovery of superconductivity in layered cobalt oxyhy- 
drate Nao.35CoO2-l.3H2O [l( has triggered intensive in- 
vestigations of superconductivity on the triangular lat- 
tice. Concerning magnetism on geometrically frustrated 
lattices, properties of antifcrromagncts with the triangle- 
based structure have been discussed for a long time [2j • In 
low-dimensional systems, the combined effects of geomet- 
rical frustration and strong quantum fluctuation cause 
peculiar behavior in low-energy physics, as typically ob- 
served in the Heisenberg zigzag chain with spin 5=1/2. 
With increasing the strength of frustration, the ground 
state is changed from a critical spin-liquid to a gapped 
dimer phase 0, 0, 0, 0| • ^ n ^ ne dimer phase, neighbor- 
ing spins form a valence bond to gain the local magnetic 
energy, while the correlation among the valence bonds is 
weakened to suppress the effect of spin frustration. 

Another important ingredient in actual materials is 
orbital degree of freedom, when electrons partially fill 
degenerate orbitals, as frequently observed in d- and /- 
electron compounds. In such a system, the interplay of 
spin and orbital degrees of freedom yields the possibility 
of orbital ordering. It is an intriguing issue to clarify the 
influence of orbital ordering on magnetic properties in 
geometrically frustrated systems. In fact, significance of 
i2 g -orbital degree of freedom has been discussed to under- 
stand the mechanism of two phase transitions in spinel 
vanadium oxides AV2O4 (A=Zn, Mg, and Cd) 0,HlIl 
It has been proposed that orbital ordering brings a spa- 
tial modulation in the spin exchange and spin frustration 
is consequently relaxed. Similarly, for MgTi2 04, the for- 
mation of a valence-bond crystal due to orbital ordering 
has been also suggested [lfl ]. 

In general, since d- and /-electron orbitals are spatially 
anisotropic, there always exist easy and hard directions 
for electron motion. Then, it is reasonable to expect that 



the effect of geometrical frustration would be reduced 
due to orbital ordering, depending on the lattice struc- 
ture and the type of orbital, in order to arrive at the spin 
structure to minimize the influence of spin frustration. 
However, the spin structure on such an orbital-ordered 
background may be fragile, since the effect of geometrical 
frustration never vanishes, unless the lattice distortion is 
explicitly taken into account. It is a highly non-trivial 
issue whether such an orbital arrangement actually de- 
scribes the low-energy physics of geometrically frustrated 
systems. In particular, it is quite important to clarify 
how the orbital-arranged background is intrinsically sta- 
bilized through the spin-orbital correlation even without 
the electron-lattice coupling. 

In this Letter, we investigate both ground- and excited- 
state properties of the e g -orbital degenerate Hubbard 
model on a ladder and a zigzag chain (see Fig. 1) by ex- 
ploiting numerical techniques. In the ladder, ferro-orbital 
ordering occurs, but the spin correlation on the rung re- 
mains finite in a paramagnetic (PM) ground state. On 
the other hand, the zigzag chain is found to be decoupled 
to a double chain in the ground state, while in the first 
spin-excited state, we find that the spin-triplet excitation 
is localized just in one side of the double chain. The ro- 
bust decoupling nature of the zigzag chain is understood 
by a concept that the multi-orbital system with geomet- 
rical frustration is spontaneously reduced to a one-orbital 
model to suppress the spin frustration at the sacrifice of 
orbital degree of freedom. This is a kind of self-organizing 
phenomenon, characteristic of the systems with a subtle 
balance among competing interactions. 

Let us consider doubly degenerate e g orbitals on the 
TV-site ladder and zigzag chain including one electron per 
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FIG. 1: Lattice location and site numbering of TV-site ladder 
and zigzag chain. The length is defined as L=N/2. 



site (quarter filling). Here the lattices are located in the 
(x,y) plane, as shown in Fig. 1. Note that the zigzag 
chain is composed of equilateral triangles. The e g -orbital 
degenerate Hubbard model is given by 

i,a,7,7',cr i,7 

+ U' ^ PiaaPiba' + J ^ AaoAbv' d ™cr> diba 
\,(t,<j' 

+J' ^ 4 7 T d i7i d i7'i d i7'T ; (!) 

where di aa (dib a ) is the annihilation operator for an 
electron with spin a in the d 3z 2_ r 2 (d x 2_ y 2) orbital at 
site i, Pi-ya—dl Jcr di lr7 , a is the vector connecting adja- 
cent sites, and t* , is the hopping amplitude between 
adjacent 7 and 7' orbitals along the a direction. The 
hopping amplitudes depend on the bond direction as 
t* a =t/A, t* b =tf a =-V3t/4, tf b =3t/4: for the x direction, 
^a=*/4, tlb=tba=V3t/4, *&=3t/4 for the y direction, 

and C=*/ 4 > *afc= < ba=v / 3^/ 8 > C= 3 */ 16 for thc oblique 
u direction. Hereafter, t is taken as the energy unit. 
In the interaction terms, U is the intraorbital Coulomb 
interaction, U' the interorbital Coulomb interaction, J 
the interorbital exchange interaction, and J' is the pair- 
hopping amplitude between different orbitals. Note that 
the relation U=U'+J+J' holds due to the rotation sym- 
metry in the orbital space and J'=J is assumed due to 
the reality of the wave function. 

In order to analyze the complex model including both 
the orbital degree of freedom and geometrical frustration, 
we employ the finite-system density matrix renormaliza- 
tion group (DMRG) method appropriate for quasi-one- 
dimensional systems with the open boundary condition 
[ill Il2| . Since one site includes two e g orbitals and the 
number of bases is 16 per site, the size of the superblock 
Hilbelt space becomes very large as m 2 xl6 2 , where m 
is the number of states kept for each block. Thus, we 
treat each orbital as an effective site to reduce the size of 
the superblock Hilbert space to m 2 x4 2 . In the present 
calculations, we keep m states up to m=200 and the trun- 
cation error is estimated to be 10~ 5 at most. Since the 
DMRG calculations have consumed CPU times, we sup- 
plementary use the Lanczos method for N=8 to study 
ground-state properties with relatively short CPU times. 

First we focus on the spin structure of the PM ground 
state at J=0, since the zigzag chain is relevant to a geo- 
metrically frustrated antiferromagnet in the spin-singlet 
PM phase. The effect of J will be discussed later. In 
Figs. 2(a) and (b), we show thc DMRG results for iV=40 
of the spin-correlation function C sp i n (i, j)=(5f Sf) with 
Sf=Y^ry(pi-yl— PiTl)/2- Note that a large value of U'=20 
is used to consider the strong-coupling region, but the 
results do not change qualitatively for small values of V. 
As shown in Fig. 2(a), we observe a simple Neel struc- 
ture in the ladder. On the other hand, in the zigzag 
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FIG. 2: The spin-correlation function measured from the 
center of the lower chain for the PM ground state in (a) the 
ladder and (b) the zigzag chain. The local magnetization for 
the first spin-excited state in (c) the ladder and (d) the zigzag 
chain. 



chain, there exists an antiferromagnetic (AFM) correla- 
tion between intra-chain sites in each of lower and up- 
per chains, while the spin correlation between inter-chain 
sites is much weak, as shown in Fig. 1(b). Namely, the 
zigzag chain seems to be decoupled to a double chain in 
terms of the spin structure. 

In order to clarify the characteristics of the spin struc- 
ture in the excited state, we investigate the local magneti- 
zation Mi=(S?) for the lowest-energy state with 5* ot =l, 
i.e., the first spin-excited state, where 5 t z ot is the z com- 
ponent of the total spin. In the ladder, the total moment 
of S'tot = l is distributed to the whole system and there is 
no significant structure, as shown in Fig. 2(c). On the 
other hand, the situation is drastically changed in the 
zigzag chain. As shown in Fig. 2(d), the total moment 
of 5f ot =l is confined in the lower chain and it forms a 
sinusoidal shape with a node, while nothing in the upper 
chain. Note that the sinusoidal shape of the local mag- 
netization is a characteristic feature of the 5=1/2 AFM 
Heisenberg chain with edges at low temperatures 0,0]. 
Thus, the double-chain nature in the spin structure re- 
mains robust even for the spin-excited state. 

Let us now consider the orbital arrangement to un- 
derstand the mechanism of the appearance of the spin 
structures. In order to determine the orbital arrange- 
ment, we usually measure orbital correlations, but due 
care should be paid to the definition. By using an an- 
gle 6\ to characterize the orbital shape at each site, we 
introduce new phase-dressed operators as |l5l | 
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FIG. 3: Optimal orbital arrangement in (a) the ladder and 
(b) the zigzag chain, (c) The angle representing the orbital 
structure in the PM ground state. The maximum value of 
T(q) at each wave vector is shown for (d) the zigzag chain 
and (e) the ladder. Note that in the zigzag chain, the wave 
vector is defined along the one-dimensional zigzag path, (f) 
The spin gap as a function of 1/L with L=N/2. 



The optimal set of {8j} is determined so as to maximize 
the Fourier transform of the orbital-correlation function, 



T(q) = (l/N^(Tfff)e 



a iq-(i-j) 



(3) 



with T?=E CT (4 )A 

In the ladder, we observe ferro-orbital (FO) ordering, 
characterized by #i/7T^1.18 in the ground state, as shown 
in Fig. 3(a). In the first spin-excited state, the FO struc- 
ture also appears, but the angle characterizing the or- 
bital shape is slightly changed as #i/7r^l.20 to further 
extend to the leg direction. On the other hand, in the 
zigzag chain, we observe that both in the ground and 
first spin-excited states, T(q) becomes maximum at q=0 
with 0j/7r~1.32, indicating a 3x 2 — r 2 orbital at each site, 
as shown in Fig. 3(b). It is emphasized here that the or- 
bital arrangement is unchanged even in the spin-excited 
state. We envisage that the orbital degree of freedom 
spontaneously becomes "dead" in low-energy states to 
suppress the effect of spin frustration. 



Consider now the effect of J on the spin-orbital struc- 
ture in the PM ground state. It is useful to show the 
results for 8-site system with the periodic boundary con- 
dition obtained by the Lanczos method in addition to a 
few DMRG results for N=40. In Fig. 3(c), we depict 
the angle 9 representing the orbital structure, where 9 is 
the optimal angle taking the same value at all sites. In 
the ladder, 9 decreases with the increase of J, indicating 
that the orbital shape is changed to further extend to 
the rung direction. Note that 9 for iV=40 is smaller than 
that for N=8. The orbital shape is sensitive to the length 
of the leg, since 9 is determined to optimize the electron 
motion both in the leg and rung directions. When we 
measure the maximum value of T(q) at each wave vec- 
tor, the amount of the FO correlation T(0,0) is reduced 
due to the effect of J, as shown in Fig. 3(d). Note that 
T(tt, 0) grows with the increase of J, while the spin cor- 
relation of q=(0, 7r) is found to develop at the same time 
(not shown here). The development of this type of the 
spin-orbital correlation indicates that the orbital degree 
of freedom is still active, since electrons occupy b orbitals 
to some extent even in the a-orbital ordered state. 

In the zigzag chain, however, no remarkable change is 
found in the angle 8 due to the effect of J, as shown in 
Fig. 3(c). In sharp contrast to the case of the ladder, the 
finite-size effect on the orbital shape is very small, since 
the orbital polarization in this case occurs due to the ge- 
ometrical reason. As shown in Fig. 3(e), irrespective of 
the value of J, T(0) takes 1/4, which is the possible max- 
imum value, indicating the perfect orbital polarization. 
All these facts suggest that the PM phase of the zigzag 
chain is described by a one- orbital model, since one of the 
doubly degenerate orbitals is selectively occupied. Thus, 
the double-chain nature of the spin structure remains ro- 
bust even for the inclusion of J. 

Let us discuss the anisotropy in the spin-exchange in- 
teractions. For the zigzag chain described by the Hub- 
bard model composed of 3x 2 — r 2 orbital, since the or- 
bital shape extends along the double chain, not along the 
zigzag path, the AFM exchange interaction along the u 
direction J\ should be much weaker than that along the 
x direction J 2 , estimated as Ji/J2=l/64 2 [l6|. Namely, 
the spin correlation on the zigzag path is reduced due 
to the spatial anisotropy of 3a; 2 — r 2 orbital. Thus, the 
zigzag chain is effectively regarded as a double-chain sys- 
tem of the 5=1/2 AFM Heisenberg chain. Concerning 
the low-energy physics, it is intuitively expected that the 
spin gap should be extremely suppressed, since the spin 
gap decreases exponentially with the increase of JW Ji in 
the gapped dimer phase in the zigzag spin chain |6j. 

On the other hand, in the ladder with the FO struc- 
ture as shown in Fig. 3(a), the orbital shape extends to 
the rung direction as well as to the leg direction. In this 
case, J rU ng/>/iog is estimated to be 0.26 0, much larger 
than Ji/J 2 = l/64 2 in the zigzag chain, where Ji cg and 
>/rung are the AFM exchange interactions along the leg 



4 




10 20 30 10 20 30 



U' IT 

FIG. 4: The ground-state phase diagram for iV=8 of (a) the 
zigzag chain and (b) the ladder. 

and rung directions, respectively. Thus, the spin corre- 
lation on the rung remains finite, leading to the simple 
Neel structure, and the spin excitation in the ladder is 
expected to be gapful similar to the spin ladder 0, 0] . 
Note that for the spin ladder of J rung / Ji eg =0.26, the spin 
gap is estimated to be 0.02 in the present energy unit. 

Let us show our numerical results on the spin gap. We 
investigate the system-size dependence of the spin gap 
A spiD =E(N/2 + 1, N/2 - l)-E(N/2, N/2) up to N=96, 
where E(N^,Ni) denotes the lowest energy in the sub- 
space with iVf up- and down-spin electrons. As shown 
in Fig. 3(f), the spin gap of the zigzag chain actually 
converges to almost zero in the thermodynamic limit, as 
expected by analogy with the zigzag spin chain. Note 
that the spin gap becomes small with the increase of J 
for each system size. On the other hand, a finite spin gap 
is observed in the ladder. Note, however, that the value 
of the spin gap is somewhat larger than that expected 
by analogy with the spin ladder, since the multi-orbital 
ladder may not be regarded as a simple spin ladder. 

Finally, let us briefly discuss the characteristics of the 
phase diagram beyond the region of small J. Here we 
show the Lanczos results for N—8, but qualitative points 
are believed to be grasped in these calculations. In 
Figs. 4(a) and (b), we show the ground-state phase dia- 
grams in the ([/', J) plane for the ladder and the zigzag 
chain, respectively. Note that the region of J>U' is ig- 
nored, since it is unphysical. Each phase is characterized 
by the total spin Stot, which is determined by compar- 
ing the lowest energies in the subspaces of S , £ )t =0~4. 
Roughly speaking, there appear three types of ground- 
state phases both in the zigzag chain and the ladder. In 
the region of small J, as discussed intensively in this pa- 
per, the ground state is the PM phase characterized by 
Stot=0. In the region of large V and J, a maximally 
spin-polarized ferromagnetic (FM) phase occurs, where 
the lowest energies in the subspaces of S* ot =0~4 are all 
degenerate. In the region next to the U'=J line, the 
ground state is characterized by iSt o t=0 again, where we 
observe a spin-spiral state, which is considered to be an 



FM state in the thermodynamic limit 20] . In addition to 
these three types of ground-state phases, we find narrow 
regions with St ot =0, 1, and/or 2. The existence of such 
a partial FM state has been discussed in a single-orbital 
Hubbard model on the zigzag chain [2l|. It is an inter- 
esting future problem to clarify the feature of the partial 
FM state in the multi-orbital system. 

In summary, we have discussed both ground- and 
excited-state properties of the e g -orbital degenerate Hub- 
bard model on the ladder and the zigzag chain. It has 
been found that the zigzag chain is reduced to a decou- 
pled double-chain spin system due to the selection of a 
specific orbital. It is considered as a general feature of 
geometrically frustrated multi-orbital systems that the 
orbital selection spontaneously occurs so as to suppress 
the effect of spin frustration. 
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